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ABSTRACT 

We consider gas at densities appropriate to protoplanetary disks and calculate 
its ability to cool due to line radiation emitted by H2O molecules within the gas. 
Our work follows that of Neufeld & Kaufman (1993; ApJ, 418, 263), expanding on 
their work in several key aspects, including use of a much expanded line database, 
an improved escape probability formulism, and the inclusion of dust grains, which 
can absorb line photons. Although the escape probabilities formally depend on a 
complicated combination of optical depth in the lines and in the dust grains, we 
show that the cooling rate including dust is well approximated by the dust-free 
cooling rate multiplied by a simple function of the dust optical depth. We apply 
the resultant cooling rate of a dust-gas mixture to the case of a solar nebula 
shock pertinent to the formation of chondrules, millimeter-sized melt droplets 
found in meteorites. Our aim is to assess whether line cooling can be neglected 
in chondrule-forming shocks or if it must be included. We find that for typical 
parameters, H2O line cooling shuts off a few minutes past the shock front; line 
photons that might otherwise escape the shocked region and cool the gas will be 
absorbed by dust grains. During the first minute or so past the shock, however, 
line photons will cool the gas at rates ~ 10^ Khr~^, dropping the temperature 
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of the gas (and most likely the chondrules within the gas) by several hundred K. 
Inclusion of H2O line cooling therefore must be included in models of chondrule 
formation by nebular shocks. 

Subject headings: ISM: lines and bands; meteors, meteoroids; shock waves; solar 
system: formation 

1. Introduction 

1.1. Radiative Transfer and Line Cooling 

In many astrophysical settings, emission of radiation via rotational and vibrational 
transitions of molecules plays an important role in cooling warm gas. Because the emission 
of such radiation is in sharp spectral lines, this mechanism is referred to as line cooling. 
Line radiation from the water molecule H2O, with its permanent electric dipole and its high 
cosmochemical abundance, is significant in a variety of settings ranging from molecular clouds 
to protostellar envelopes (e.g., Cernicharo & Crovisier 2005). More recently, line cooling 
from H2O molecules has been recognized to play a pivotal role in the energetics following 
the passage of a shock wave through the dense gas in the solar nebula protoplanetary disk 
(see Desch et al. 2005). 

The rate of emission of line radiation from a warm gas containing water molecules is 
difficult to calculate, because there are so many accessible rotational and vibrational energy 
levels, and therefore a great many transitions. In dense molecular gas (nH2 > lO^'^cm"'^), 
these energy levels are populated according to Boltzmann statistics (i.e., in local thermody- 
namic equilibrium, or LTE), but at lower densities the populations must be calculated by 
balancing transition rates. Finally, the cooling of gas by line radiation relies on the ability of 
the photons to escape the system without being reabsorbed; to the extent that other, nearby 
molecules can absorb the emitted photons, the gas does not cool. Photons generated during 
a molecular transition are emitted and also reabsorbed over a small range of frequencies 
centered on the line frequency, with an efficiency that depends on the frequency shift from 
the line center. Rather than integrating the equations of radiative transfer over all frequen- 
cies within the line, all standard treatments of the problem instead assume integration over 
this "line profile" and use a frequency-integrated escape probability of photons. This escape 
probability is a function of the column density of the molecule within the system. This is 
the approach taken by Neufeld & Kaufman (1993; hereafter NK93), in particular. 

The calculation of NK93 of the rate of line cooling from H2O molecules has remained 
the state of the art for 15 years; but, as we discuss below, several aspects of the NK93 
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calculation are now out of date. Their calculation uses a limited database of transitions, 
an oversimplified escape probability formulism, and ignores absorption by dust grains. It is 
one of the goals of the work here to update the H2O line cooling rates to improve on the 
calculations of NK93. A second goal of this paper is assess whether or not line cooling plays 
a significant role in cooling the gas following the passage of a shock through the dense gas of 
the solar nebula. As discussed by Desch et al. (2005), previous modeling has not determined 
whether or not line cooling can be neglected. This has significant implications for the issue 
of chondrule formation. 



1.2. Chondrules 

The parent bodies of the most primitive meteorites, the chondrites, formed at about 
2-3 astronomical units (AU) from the Sun, 4.57 billion years ago (Wadhwa & Russell 2000). 
Chondrites are the most primitive meteorites in our collection, in that they have suffered very 
little alteration since their formation, and therefore contain information about the conditions 
that existed in the early solar nebula. Chondrites are remarkable for containing calcium- 
rich, aluminum-rich inclusions (CAIs), the oldest solids in the Solar System, whose formation 
has been dated to 4567.2 ± 0.6 Mya (Amelin et al. 2002). Also found in abundance within 
chondrites are sub-millimeter- to millimeter-sized, (mostly ferromagnesian) igneous spheres, 
called chondrules, from which the chondrites derive their name. Chondrules formed, at 
most, ~ 2-3 million years after CAIs (Amelin et al. 2002; Kita et al. 2005), as melt droplets 
that were heated to high temperatures while they were independent, free-floating objects 
in the early solar nebula. After they were heated, cooled and recrystallized, chondrules 
were incorporated into the parent bodies from which chondrites originate. Chondrules are 
capable of providing incredibly detailed information about conditions in the Solar System 
protoplanetary disk, if the process that led to their heating, melting and recrystallization 
could be understood. Chondrules make up to 80% of the volume of ordinary chondrites 
(Grossman 1988), and it is estimated that ~ 10^^ g of chondrules exist in the asteroid belt 
today (Levy 1988). Such a prevalence of chondrules suggests that chondrule- forming events 
were widespread in the solar nebula. A process that can melt 10^^ g of rock is surely a 
dominant process in the solar nebula disk. 

The chondrule formation process, despite its obvious importance, has been a mystery 
in the field of meteoritics for two centuries (Sorby 1877). Proposed mechanisms include 
interaction with the early active Sun, through jets (Liffman & Brown 1995; Liffman & 
Brown 1996) or solar flares (Shu et al. 1996, 1997, 2001), melting by lightning (Pilipp et al. 
1998; Desch & Cuzzi 2000), and melting by planetesimal impacts (Urey & Craig 1953; Urey 
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1967; Sanders 1996; Lugmair & Shukolyukov 2001). The dominant theory, though, is that 
chondrules were melted in shock waves in the protoplanetary disk (Hewins 1997; Jones et al. 
2000; Connolly & Desch 2004; Connolly et al. 2006). A shock wave, or shock front, is a sharp 
discontinuity between supersonic and subsonic gas, over an area only a few molecular mean 
free paths thick, typically only meters in the solar nebula. The gas is slowed, compressed, 
and heated by the time it reaches the other side of the shock front. Solids moving with the 
gas are heated not only by thermal exchange upon entering the shocked region, but also 
by friction as they are slowed to the post-shock gas speed, and by absorbing the infrared 
radiation emitted by other solids. The shock model of chondrule formation appears able 
to resolve the chondrule formation mystery, because it makes several detailed predictions 
about chondrule formation that are largely borne out by observation and experimentation, 
especially regarding chondrule thermal histories. 



1.3. Chondrule Formation by Nebular Shocks 

Heating and cooling rates of chondrules have been determined experimentally. Accord- 
ing to furnace experiments, in which melt droplets with chondrule compositions are allowed 
to cool and crystallize, reproduction of chondrule textures requires particular ranges of cool- 
ing rates between the liquidus temperature (~ 1800 K) and solidus temperature (~ 1400 
K; Hewins & Connolly 1996). Porphyritic chondrules cooled at about 10 — lO^Khr"^, and 
barred-olivine chondrules cooled at about lO^Khr"^ (Hewins et al. 2005; see also Desch & 
Connolly 2002). Additionally, chondrules retain volatile elements such as S, indicating that 
they did not remain above the liquidus for more than minutes, and cooled quite rapidly 
(> lO^Khr"^ while above the liquidus (Yu & Hewins 1998)). Finally, there is no indication 
of the isotopic fractionation that would arise from the free evaporation of alkalis such as Na, 
which constrains the time spent at high temperature before melting (Tachibana et al. 2004). 
Modeling of isotopic fractionation has shown that chondrules must heat up from 1300 to 
1600 K on the order of minutes or less (Tachibana & Huss 2005). 

Passage through nebular shocks satisfies nearly all the experimental constraints on chon- 
drule formation in broad brush (lida et al. 2001, hereafter INSN; Desch & Connolly 2002, 
hereafter DC02; Ciesla & Hood 2002, hereafter CH02; Desch et al. 2005). Prior to passage 
through the shock front, chondrules are moderately heated by absorbing radiation emitted 
by chondrules which have already passed through the shock front. Upon passage through the 
shock front, the gas is immediately slowed, and compressed and heated, but the chondrules 
continue at supersonic speeds through the gas. They achieve their peak heating during this 
stage, due to absorption of radiation, thermal exchange with the gas, and supersonic fric- 
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tional drag heating. This stage lasts until the chondrules slow to the gas velocity, which 
takes an aerodynamic stopping time, tgtop = PsO-s/PgC ~ 1 minute (where ps and are the 
particle density and radius, and Pg and C are the post-shock gas density and sound speed). 
At this time, the gas and chondrules are dynamically coupled and chondrules are heated only 
by absorption of radiation and thermal exchange with the gas. Soon thereafter the solids and 
gas become thermally coupled as well, and both components achieve similar temperatures. 
Although the chondrules pass through the first 100 km of the post-shock region more rapidly 
than the gas, once they do, they will thermally equilibrate to the gas temperature. Once this 
happens, the gas and chondrules cool together as fast as they can either radiate away energy 
or leave the source of infrared radiation that heats them. DC02 showed that the cooling rate 
of gas and chondrules should be 

CR = 50 ( ^ -] ( (s+-] KhT~\ (1) 

Vl0-9gcm-V V7kms-V V 50/ ' 

where pg is the pre-shock gas density, Vs the shock speed, C is the "concentration" of 
chondrules (the chondrules-to-gas mass ratio normalized to 3.75 x 10^^, where a chondrule 
radius 300 pm has been assumed) , and 6 the concentration of submicron dust (the dust-to- 
gas mass ratio normalized to 1.25 x 10""^). More precisely, S should be interpreted as the 
Rosseland mean opacity of the dusty gas, normalized to 1.14 cm^ g~^. For typical parameters, 
then, cooling rates ~ 10^ Khr~^ are obtained, consistent with the cooling rates of chondrules. 



1.4. Estimates of Line Cooling in Shocks 

The effect of line cooling — the cooling of gas by emission of so-called line photons by trace 
molecules in the gas such as CO and H2O — in solar nebula shocks has been considered by 
INSN and by Miura & Nakamoto (2006). In both cases the cooling rates of NK93 were used. 
INSN assumed a gas optically thin to the line radiation and found chondrule cooling rates 
~ lO^Khr"^ in many cases. Miura & Nakamoto (2006) allowed the gas to become optically 
thick to this line radiation, using the formulation of NK93 described below. They found that 
for plausible parameters the cooling rates of chondrules and gas cluster around 5000 Khr~^, 
with line cooling playing an important role in the cooling of the gas. Chondrules and gas are 
thermally coupled only a few minutes after passing through the shock front (INSN; DC02; 
CH02), so these chondrule cooling rates are controlled by the rate at which gas cools by 
emission of line radiation. To understand the essence of their results without reproducing 
all of their calculations in detail, we now estimate the cooling rate of chondrules and gas 
following the shock, using the tabulations of NK93. 

We first examine how the cooling rates predicted by NK93 and used by INSN and Miura 
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& Nakamoto (2006) are calculated. For a wide range of densities and temperatures, NK93 
calculated the steady-state level populations of H2O and CO molecules (denoted 'M'), and 
the rates of emission of photons. The escape of these photons from a system was considered 
under the large-velocity-gradient (Sobolev) approximation as a function of column density 
n{M)d, via use of the parameter 

^ - Til- 

These results were adapted to a large range of geometries of interest, including those not 
consistent with the Sobolev approximation, in particular the one of relevance to nebular 
shocks, a static, plane-parallel slab of thickness d. According to NK93, the cooling rate at 
the center of the slab (at depth d/2) will depend on the similar parameter 

N^'<^, (3) 
Av 

where Av is the thermal velocity of the molecules M. An important assumption made by 
NK93 is that if the line photons can escape to the edge of the slab, then the gas is cooled; 
if the photons do not escape to the edge of the slab, they are essentially re-absorbed on 
the spot and there is no cooling. The cooling rate per volume is then n{m)CuTE{N), where 
the subscript 'LTE' is the notation used by NK93 to denote cooling rates derived under the 
assumption that the population levels follow a Boltzmann distribution, which is the regime 
of interest for solar nebula shocks. NK93 have tabulated £lte(A^) for both vibrational and 
rotational lines of H2O and CO, for every decade in A^. 

The plane-parallel slab is equivalent to the nebula shock geometry because of the large 
velocity jump at the shock front. In a typical chondrule-forming shock, the gas velocity 
changes nearly instantaneously from ^ 7kms~^ to ^ 1.6kms~^, a jump of over 5kms~^. 
This velocity difference is much larger than the thermal velocities of the H2O molecules 
[Av = {2kT/mn2of''^ = 1.36 km s"^ at 2000 K] that control the width of the line profile. 
Essentially, once a line photon emitted from the post-shock gas escapes to the shock front, it 
will continue to travel relatively unimpeded by gas. Only one change to the NK93 formulation 
must be made before it is adopted to the nebular shock geometry. For gas a distance z past 
the shock front, the equivalent thickness of the slab is li = 2z, and the cooling rate is only 
half the cooling rate calculated by NK93, who assumed radiation could escape either side of 
the slab. 

The cooling rate of the gas yields 

^ ^ 2.8nH2A:?^ = -^r2(H20)£LTE(iV), (4) 



dt dt \-f-lJ ' dt 2 

where we have ignored compression of the gas after passing through the shock front in this 
toy model, and we have assumed an abundance of He of 10% by number. Converting the 
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time derivative to a spatial derivative by assuming a constant gas velocity Vg past the shock 
front, we find 

dT _ Av £LTE(iV) 

dN b.QrimVgk 2 ' ^ ' 

where we have assumed a constant Av throughout the post-shock region. In actuality, Av 
would decrease past the shock front, but that would merely reduce the total cooling, so our 
assumption overestimates the degree of cooling. This formula provides the gas temperature 
as a function of A^, the effective column density past the shock front. 

Fedkin & Grossman (2006) showed that the rimo/nm ratio in the solar nebula probably 
varied with temperature, but found that ?t,h2o/''^H2 = 5 x 10~^ was a good average for the 
temperatures of interest. We have adopted a canonical ratio ?t-h2o/'^H2 = 8 x 10~^. This 
is only slightly higher than some commonly assumed canonical ratios but is consistent with 
the abundance of Lodders (2003), nH2o/''^H2 = 8.88 x 10~^. We note that the exact value of 
nH2o/'^H2 in the chondrule formation environment is likely variable with time (Ciesla & Cuzzi 
2006), and consider variations in this ratio later in the paper. Some dissociation of H2O is 
likely to occur at the peak post-shock temperatures we consider, but we are interested in the 
maximum cooling possible by line photons, so we neglect dissociation in our calculations. 

Interpolating between the tabulated values of £lte(^) provided by NK93, we have 
integrated equation [S] to find the total drop in temperature by the time the gas is at an 
equivalent of iV = 10^^ cm~^ km~^ s past the shock front. For typical parameters (pre- 
shock density rim = 2 x lO^^cm"^ and velocity Vg = 7kms~^; mixing ratio nmo/^m = 
8 X 10~^; post-shock density rim = 1-2 x lO^^cm"'^ and velocity Vg = 1.2kms~^; and 
post-shock thermal velocity Av = 1.6 km s"^), this equates to d = 1.3 x 10^ km, a distance 
z = 6.6 X 10'^ km, and a time 1.6 hours after passing through the shock front. (Again, this toy 
model neglects the compression and slowing of the gas. It also assumes wrongly that Av is 
constant, whereas it should decrease slightly as the gas cools. It is nonetheless illustrative.) 
At these distances past the shock front, the total cooling due to rotational lines of H2O 
and CO is completely negligible, < 4K. This illustrates that the rotational lines become 
optically thick so quickly that they are unable to escape the post-shock region and cool 
the gas. On the other hand, the total cooling from CO vibrational line photons is roughly 
72 K. This leads to an average cooling rate of ~ 40Khr~^, which is important, but not 
significantly greater than the cooling rates of chondrules (~ lO^Khr"^). The assumption 
of INSN that the lines of H2O and CO are optically thin is therefore invalidated, because 
these lines are indeed optically thick. The final line cooling mode (not considered by INSN, 
but considered by Miura & Nakamoto 2006) is vibrational line cooling from H2O. We find a 
total cooling due to these line photons on the order of 800 K. This is a significant cooling of 
the gas: the average cooling over 1.6 hours is ~ 500Khr-\ with even higher cooling rates 
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obtained early on. The cooling of the gas by vibrational line photons of H2O, despite being 
somewhat optically thick, is the source of the high chondrule cooling rates predicted by 
Miura & Nakamoto (2006). We confirm the statement made by Miura & Nakamoto (2006) 
that "...line emission is important for the gas cooling." We hereafter focus on cooling due 
to H2O alone, as the dominant coolant is expected to be water. 

Clearly, any model of chondrule cooling rates in solar nebula shocks must account for 
the emission of line photons from H2O. However, incorporating the cooling rates of NK93 
directly into a shock code, as Miura & Nakamoto (2006) did, is not ideal. First, line cooling 
due to rotational and vibrational transitions of H2O was calculated by NK93 using roughly 
50,000 transitions from the HITRAN database (Rothman et al. 1987). While this was the 
best available data at the time, much more extensive databases now exist. Second, NK93 
calculated escape probabilities under the large-velocity-gradient (Sobolev) approximation, 
and asserted that the results will apply to the case of a static, plane-parallel slab if is 
defined as described above. Third, NK93 assumed a one-sided escape probability for line 
photons equal to 0.5/(1 -|- 3r), where r is the Sobolev optical depth; more exact escape 
probabilities exist, as described below. Fourth, it is not necessarily the case that if photons 
do not escape to the edge of the slab that the gas does not cool locally; photons re-absorbed 
halfway to the edge still cool the gas at the slab center. Finally, NK93 did not consider the 
absorption of line photons by intervening dust. Inclusion of dust is absolutely mandatory, as 
the following argument makes clear: For a solar-composition gas with 0.5% of the mass in the 
form of a = 0.5 fim radius grains (similar to the size of matrix dust grains in meteorites), the 
opacity to short- wavelength radiation is as high as k = (pdust/Pgas) (3/4psa) = 30cm^g~^, 
where the internal density of the dust grains is ps ~ 2.5 gcm"'^ . Over the distance for which 
the optical depth in dust is unity, the parameter N is therefore less than 10^^ cm~^ km~^ s, 
assuming a ratio 7T'H2o/''^H2 = 8 x 10~^. This corresponds to distances ~ 10"^ km, or times of 
about 10 minutes, after passage through the shock front. As this is smaller than the A^ for 
which the most significant cooling takes place, dust grains are capable of absorbing photons 
that would otherwise cool the gas. A thorough treatment of line cooling should improve on 
the calculations of NK93 in the areas identified above. 



1.5. Outline 

In this paper we calculate the cooling rate of gas due to emission of line photons by 
H2O molecules, improving upon the work of NK93 with the use of an expanded database, 
improved escape probabilities, and the inclusion of dust. First we evaluate escape probability 
approximations, both in the presence and absence of dust. We then calculate the cooling 
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rates from H2O molecules using the 1.2-million line SCAN-H2O database (Jorgensen et al. 
2001). Finally, we incorporate these cooling rates into a toy model of chondrule thermal 
evolution. We conclude that for canonical parameters, line cooling is significant in the first 
few minutes past the shock front, but that afterward dust grains absorb the line photons 
that would otherwise cool the gas. We discuss the effect of the important dust-to-water ratio 
on these conclusions. 



2. Calculation of H2O Cooling Rates 

We have used the SCAN-H2O database of J0rgensen et al. (2001) of 1.2 million lines 
(rotational plus vibrational) of the H2O molecule to calculate the cooling rates due to emis- 
sion of line photons from H2O. We continue to assume, as NK93 did, that if photons are 
not able to escape the region entirely, they are re- absorbed on the spot. This assumption 
will be relaxed in future work, but for now it allows a crude estimate of the importance of 
line radiation. In other respects, we improve on the NK93 calculation by more accurately 
calculating escape probabilities, and by including the possibility of absorption by dust grains. 



2.1. Escape Probabilities Without Dust 

The probability that a line photon will escape a semi-infinite volume, otherwise known 
as the one-sided escape probability, Pesc, is given by 

Pesc = ^ / ^X) E2{t^x)) dx, (6) 

^ A-Ao 

AA^ ' ^ ' 



c 2kT 
AXd = — 



Ao V fT^mo 

(Avrett & Hummer 1965; Hummer & Rybicki 1971; Rybicki & Lightman 1979; Bowers & 
Deeming 1984). Here $(x) is the line profile, x is the frequency measured from line center 
in Doppler widths (AAd), and E2 is the second exponential integral. In equation [6l r refers 
to the optical depth integrated over the line, which is a/tt times the optical depth at line 
center. Considering the combined effect of both Doppler and Lorentz broadening, the line 
profile is given by the normalized Voigt profile 

*(^) = ^ / 7 ^ (9) 
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(Avrett & Hummer 1965; Rybicki & Lightman 1979; Bowers & Deeming 1984), where a is 
the ratio of Lorentz to Doppler width (Avrett & Hummer 1965). In astrophysical situations 
a ^ 1, perhaps as large as 0.1 (Mihalas 1978; Bowers & Deeming 1984). In this study, the 
Voigt profile was calculated using the algorithm of Zaghloul (2007), where the Voigt function 
is written as a single proper integral with a damped sine integrand. A straightforward 
numerical integration was performed to calculate the "exact" one-sided escape probability 
Pesc in the absence of dust. 



Besides the form adopted by NK93, several other approximations to the one-sided es- 
cape probability exist in the literature. Hollenbach & McKee (1979) give the following 
approximation to the one-sided escape probabihty: 

^csc = - ^ 77:7 • (10) 

2 l + r(2 1n(2.13 + r77r))^/' 
The approximation given by Collin- Souffrin et al. (1981) is 

_ 1 1 

-'esc 1 /o • (11) 

2 l + 2r(ln(r/0F+l))'/' 



Mathews (1992) estimates 



esc „ 1 /o : 

2 l + 2r(ln(rv/?+l))'/' 

and Dumont et al. (2003) give 



(12) 



1 - 

Pesc = — , r < 1, (13) 

4r 

^ , r>l. (14) 



2v/?r(l.2 + ln^f5^ 

These approximations are all in the Doppler limit, in which a = and the line profile is 
given by 7r~^/^e~^^. We have evaluated the approximations to the escape probabilities of 
Hollenbach & McKee (1979), CoUin-Souffrin et al. (1981), Mathews (1992), NK93, and 
Dumont et al. (2003), and compared the results to our "exact" escape probabilities (Figure 
1). The maximum error between these approximations and the exact escape probability 
are as follows: 11.42% for Hollenbach & McKee (1979); 15.70% for Collin- Souffrin et al. 
(1981); 31.34% for Mathews (1992); and 37.22% for Dumont et al. (2003). The maximum 
error between the approximation used by NK93 and the exact escape probability is 46.44%, 
although it is not quite appropriate to make this comparison, as NK93 used the Sobolov 
optical depth in their calculations. The approximation of Hollenbach & McKee (1979) is 
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Fig. 1. — Escape probability approximations used by Hollenbach & McKee (1979), CoUin- 
Souffrin et al. (1981), Mathews (1992), NK93, and Dumont et al. (2003), compared to the 
"exact" escape probabihty calculated here (see text). 
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superior to the others and provides a satisfactory fit; we will use this approximation in what 
follows. 

Using the HoUenbach & McKee (1979) approximation to the escape probability, a param- 
eter study was conducted to determine the effect of varying the value of a in the calculation 
of the Voigt function. It was found that varying the value of a between 0.0 and 0.1 increased 
the escape probability significantly for r ^ 1 (Figure 2), consistent with figures 4 and 5 
of Hummer & Rybicki (1982). It is not immediately clear how this will affect the overall 
cooling rate, but we anticipate a small effect as the lines at which the escape probability 
increases with increasing a are optically thick, leading to lower cooling rates. We investigate 
this effect of the Voigt function on the cooling rates in Section 3. 



Dust can also absorb photons, preventing their escape. Hence, the proper treatment of 
the escape probability must account for the absorption by dust, where the optical depth to 
line photons is now given by + r$(x), and the exact escape probability is 



This equation reflects the fact that dust grains provide a continuum opacity that is always 
capable of absorbing photons. We have used this equation to calculate the exact escape 
probabilities in the presence of dust (Figure [3]). Because it is computationally burdensome 
to calculate the escape probability for all possible combinations of r and Td, we seek an 
approximation such that the escape probability is the product of two separate, independent 
functions of r and t^. On physical grounds, one expects this function to resemble Pcsdj, Td) ~ 
-Pesc(T"d = 0) X _E'2(rd). Clearly this is the appropriate expression when Td ~ 0. It is also exact 
in the limit r = 0. The deviation is expected to be strongest when r ~ Td ~ 1, but happily 
the deviation is not great. Applying the slight modification. 



provides an excellent fit to the exact probability calculated using equation [151 as illustrated 
in Figure [3], and tabulated in Table 1. 

It can be seen that the approximations including dust fit very well in all regions except 
where r ~ 1 and Td ~ 1, as expected. We can apply these escape probabilities to the 
calculation of the cooling rate from water when dust is present. 



2.2. Escape Probabilities With Dust 




(15) 



n:r(^,^d) = Pcsc(rd = 0) X E2(0.95rd), 



(16) 
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Fig. 2. — Exact escape probability for three different values of a (the ratio of Lorentz to 
Doppler width). 

Table 1. Maximum Error for Pcsc, with the inclusion of dust. 





maximum error (%) 


0.0 


11.42 


0.03 


24.52 


0.1 


25.02 


0.3 


47.51 


1.0 


54.27 


3.0 


47.73 


10.0 


65.20'^ 


30.0 


128.04^^ 



^Maximum error was 
reached at high optical depth 
where the escape probability 

< 10-8. 
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Fig. 3. — Approximations to the escape probability with the inclusion of dust (dashed lines), 
compared to exact escape probabilities, including dust (solid lines). From left to right, the 
values used for t^/tq are 30, 10, 3, 1, 0.3, 0.1, 0.03, and 0. 
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2.3. Line Cooling 

The cooling per line in the optically thin limit is given by 

^ni _c.m SnkT f hc/XkT\ 

where S{T) is the temperature-dependent line strength of the spectral line, and A is the 
wavelength at line center (see Appendix). In our calculations, we assume that level popula- 
tions obey a Boltzmann distribution. This is justified if the number density of protoplanetary 
disks greatly exceeds the critical density of each transition. As NK93 found no variations in 
level populations when the density exceeded ~ 10^^° cm~^, we assume that all level popula- 
tions have critical densities below this value. Since we are motivated by chondrule formation, 
for which the gas densities exceed ~ 10^^ cm""^, the use of a Boltzmann distribution is jus- 
tified. Unfortunately, for other problems (e.g., molecular clouds), our approach is limited in 
applicability. 

The cooling rate is reduced below its optically thin limit as the column densities are 
increased, due to the inability of photons to escape the system. Specifically, the total cooling. 
Cute, is given by the total cooling per line, summed over all lines, including the escape 
probabilities: 

ClTE = Kl PcsciTuh 1"d)- (18) 

Based on the methods of Plume et al. (2004), the gas optical depth at line center, tq, can 
be found from the linestrength and column density of water: 

[Av/c) V 

(see Appendix), where the Doppler hnewidth is given by 

( 1kT\^ 

Av = , (20) 



and is ~ 1.36 km s~^ at T = 2000 K. In terms of the optical depth at line center, r^i = TqX y^. 
The dust optical depth at a depth z into a semi-infinite volume is given by 

Td=pgZK,{X), (21) 

where k(A) is the opacity per gram of gas. Assuming a dust-to-gas ratio Pd/Pg = 5 x 10^^, 
a particle radius = 0.5 pm, and an absorptivity Qabs = 1 for A < 2TTas and Qabs = Inas/X 
for A > 27ras, we derive 

K = 30 min [l, (A/3.1 pm)'^] cm^ g-^ (22) 
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We have determined the Rosseland mean dust opacity to be = 28.76 cni^g"^ at T = 
2000 K. Our estimates are similar to those derived by Henning & Stognienko (1996) using a 
particle size distribution instead of our simplified monodispersion. 

2.4. Cooling Rates 

We first calculate the cooling rates in the absence of dust. Using many more transitions 
(and improved escape probabilities) we find that, in the absence of dust, our cooling rates are 
enhanced by about 30% over those calculated by NK93 (Figure H]). For the cooling without 
dust, we found a slight difference (< 6%) in cooling rate between the cases where the escape 
probabilities are calculated exactly, and when we use the approximation of HoUenbach & 
McKee (1979). For consistency, we will quote cooling rates obtained using the HoUenbach 
& McKee (1979) approximation. The cooling rate for the T = 2000 K case, in the absence 
of dust, is plotted as the rightmost curve in Figure O Also plotted in Figure \5\ are the 
cases where the dust-to-water ratio is varied from 300 times its canonical value to 1/300 
times its canonical value. Here the canonical case refers to a dust-to-gas mass ratio of 0.5%, 
yielding the opacity discussed above; the dust-to-water ratio effectively measures only the 
dust opacity per gram of water vapor. The canonical amount of water assumed here is a 
ratio H2O/H2 = 8 X 10~^. A higher dust-to-water ratio means dust grains are more likely to 
absorb line photons, and the cooling rate due to line emission is reduced below the optically 
thin limit at a smaller total column density of water. Surprisingly, even small amounts 
of dust (300 times smaller than the canonical limit, or a dust-to-gas ratio of ~ 0.001%) 
will prevent the gas from cooling as freely as it would without dust at column densities 
A'h20 > lO^^cm"^. The maximum error in each case between the exact value of the cooling 
rate (using the escape probability of equation 15) and that found using the approximation 
to the escape probability (equation 16) is shown in Table 2 (for T = 2000 K). Figure [6] gives 
the same information for the case when T = 1500 K. The cooling rates (assuming a = 0, the 
pure Doppler broadening case) for T = 1250 K, 1500 K, 1750 K, 2000 K, and 2250 K, both 
in the absence of dust and with the canonical abundance of dust, are tabulated in Tables 3 
through 12. 

Absorption of line photons by dust grains drastically affects the cooling of the gas. The 
dust opacity is a strong function of wavelength, so we investigate the effects of assuming 
a constant opacity. In the continuum opacity case it would be reasonable to assume the 
flux of photons propagating through a gas is equivalent to the flux derived assuming the gas 
had a wavelength-independent opacity equal to the Rosseland mean opacity. We therefore 
recalculated the cooling rate when the dust opacity equaled the Rosseland mean opacity 
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Table 2. Maximum Discrepancy between "Exact" and Approximate Cooling Rates shown 

in Figure El at T = 2000 K. 



Dust-to- Water Ratio*^ 


maximum error (%) 


0.0 


5.54 


1.0 


14.01 


3.0 


13.29 


10.0 


12.24 


30.0 


11.04 


100.0 


9.26 


300.0 


7.45 


1/3 


14.77 


1/10 


15.64 


1/30 


15.39 


1/100 


17.93 


1/300 


18.40 



Relative to the canonical value. 
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cited above, at all wavelengths. For the canonical dust-to-water ratio, the results are plotted 
in Figure [71 The optical depth of dust is given by 

Td = 1.4mH2 (23) 

nmo/nm 

Using the Rosseland mean opacity and canonical water-to-gas ratios cited above, we would 
estimate r ~ 1 when A^h20 ~ 6 x 10^^ cm~^. This is indeed the point where the cooling 
rate is significantly reduced due to absorption by dust grains. When the opacity is replaced 
by the wavelength-dependent opacity however, the cooling is not effectively reduced 

until higher column densities are reached. This signifies that most of the cooling is effected 
by emission of line photons with wavelengths at which /t(A) < kr, i.e., A > 3 fim. Use of a 
wavelength- dependent opacity in conjunction with a calculation of the cooling at each H2O 
wavelength is therefore necessary. We use the wavelength-dependent opacity in all cases 
cited here. 

Finally, we investigate the effect of the line broadening parameter a on the cooling 
rate. Increasing a above the a = pure Doppler broadening case has the effect of putting 
more emission in the optically thin wings of the line profile. Effectively, this should act 
like a reduction of the overall water column density. We have calculated the cooling rates 
assuming line profiles with a = and a = 0.1, in the cases where no dust is present, and 
when it is present at the canonical value. These cooling rates are plotted in Figure [HI As 
expected, the cooling rate is slightly higher when a = 0.1, because the column density of 
water is effectively reduced; however, the difference is slight. In the case with no dust, the 
discrepancy between the a = and a = 0.1 cases is < 11% for a column density of water of 
Nmo < lO^^cm"^, increasing to ~ 41% at a column density of lO^^cm"^. In the case with 
dust, the discrepancy is < 13% at A^h20 < lO^^cm"^. These calculations are much more 
computationally burdensome than the a = pure Doppler broadening case, because of the 
need to calculate the Voigt function. Fortunately, it is seen that their effects on the cooling 
rates of gas are minimal, < 10% in the regime where cooling is significant, for a up to 0.1. 
As a is typically ^ 1 (Mihalas 1978; Bowers & Deeming 1984), we will assume a = in the 
cases that follow. 



3. Cooling During Chondrule-forming Shocks 

Line cooling is effective at cooling gas, provided the line photons are not reabsorbed 
by H2O molecules or dust grains before they can escape to the cool, pre-shock region. In a 
nebular shock, it is assured that sufficiently far from the shock front, no line photons can 
escape and the gas and dust will become thermally coupled and cool slowly, if at all. What 
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Fig. 4. — The cooling rates of NK93 compared to those of this study. 

Table 3. Cooling rates without the inclusion of dust, a = 0.0, T = 1250 K 



H20 



C 



LTE exact 



c 



LTEa 



70 error 



1013 


1 


547 


X 


10- 


-12 


1.548 


X 


10- 


-12 


0.03 


1014 


1 


543 


X 


10- 


-12 


1.546 


X 


10- 


-12 


0.16 


1015 


1 


517 


X 


10- 


-12 


1.526 


X 


10- 


-12 


0.61 


1016 


1 


405 


X 


10- 


-12 


1.425 


X 


10- 


-12 


1.44 


lOi^ 


1 


128 


X 


10- 


-12 


1.160 


X 


10" 


-12 


2.82 


1018 


6 


259 


X 


10" 


-13 


6.415 


X 


10" 


-13 


2.95 


1019 


2 


227 


X 


10- 


-13 


2.293 


X 


10" 


-13 


2.96 


1020 


6 


611 


X 


10" 


-14 


6.816 


X 


10- 


-14 


3.10 


1021 


1 


615 


X 


10" 


-14 


1.668 


X 


10" 


-14 


3.33 


1022 


3 


318 


X 


10" 


-14 


3.437 


X 


10" 


-14 


3.58 


1023 


5 


750 


X 


10- 


-16 


5.978 


X 


10- 


-16 


3.95 


1024 


8 


587 


X 


10- 


-17 


8.959 


X 


10- 


-17 


4.34 
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Fig. 5. — Exact cooling rates due to H2O for various dust-to-water ratios (solid lines), as a 
function of water column density, and the cooling rates calculated using the approximation 
to the escape probability (dashed lines) for T = 2000 K. The rightmost curve is the case 
with no dust; the red curve is the canonical case (see text). The curves to the left of the 
canonical curve (from right to left) show cooling rates for 3, 10, 30, 100, and 300 times the 
canonical dust-to-water ratio. The curves to the right of the canonical curve (from left to 
right) are 1/3, 1/10, 1/30, 1/100, and 1/300 times the canonical dust-to-water ratio. 
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Fig. 6. — Same as Figure [5], but with T = 1500 K. 

Table 4. Cooling rates with the canonical dust-to-water ratio (see text), a = 0.0, 

T = 1250 K 



Nh20 






•^LTEapprox 




% error 


10^3 


1.547 X 10" 


-12 


1.548 X 10" 


12 


0.03 




1.543 X 10- 


-12 


1.546 X 10" 


12 


0.16 




1.516 X 10" 


-12 


1.525 X 10- 


12 


0.59 


10^6 


1.397 X 10" 


-12 


1.415 X 10" 


12 


1.26 


10^^ 


1.083 X 10' 


-12 


1.097 X 10" 


12 


1.62 




4.686 X 10- 


-13 


4.439 X 10- 


13 


5.27 




3.204 X 10- 


-14 


2.765 X 10- 


14 


13.71 


1020 


1.459 X 10- 


-17 


1.226 X 10- 


17 


15.98 


1021 


1.101 X 10- 


-28 


2.148 X 10- 


28 


95.04 


1022 


0.00 




0.00 




0.00 


1023 


0.00 




0.00 




0.00 


10^4 


0.00 




0.00 




0.00 
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Fig. 7. — Cooling rates due to H2O using the Rosseland mean dust opacity at T = 2000 K 
(dashed hne) and using the wavelength-dependent dust opacity (solid line). Most line cooling 
takes place at wavelengths longward of 3 /im, where the dust is relatively optically thin. Use 
of a single dust opacity at all wavelengths is not warranted. 
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Fig. 8. — Cooling rates with a broadened line profile corresponding to a = 0.1 (solid curves), 
compared to those calculated assuming pure Doppler broadening a = (dashed curves), 
both in the absence of dust (top curves) and with a canonical dust-to-water ratio (bottom 
curves) . 
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Table 5. Cooling rates without the inclusion of dust, a = 0.0, T = 1500 K 



Nh20 -^LTEexact -^LTEapprox ^ ^1101 



1 ni3 

iU 


9 

Z, 


7n/i 


X 


iU 


-12 


o 

Z, 


. ( uo 


X 


lU 


-12 


n 


no 
.uz 




2, 


.700 


X 


10" 


-12 


2, 


.704 


X 


10- 


-12 


0, 


.12 


10^5 


2, 


.667 


X 


10" 


-12 


2, 


.679 


X 


10- 


-12 


0, 


.47 




2, 


.519 


X 


10" 


-12 


2, 


.548 


X 


10- 


-12 


1, 


.15 


10^^ 


2, 


.127 


X 


10- 


-12 


2, 


.179 


X 


10- 


-12 


2, 


.47 




1, 


.347 


X 


10- 


-12 


1, 


.383 


X 


10- 


-12 


2, 


.81 




5, 


.597 


X 


10- 


-13 


5, 


.755 


X 


10- 


-13 


2, 


.84 


1020 


1, 


.850 


X 


10- 


-13 


1, 


.905 


X 


10- 


-13 


3, 


.00 


1021 


4, 


.992 


X 


10- 


-14 


5, 


.153 


X 


10- 


-14 


3, 


.26 


1022 


1, 


.094 


X 


10- 


-14 


1, 


.132 


X 


10- 


-14 


3, 


.45 


1023 


1, 


.917 


X 


10- 


-15 


1, 


.991 


X 


10- 


-15 


3, 


.86 


1024 


2, 


.715 


X 


10- 


-16 


2, 


.837 


X 


10- 


-16 


4, 


.48 



Table 6. Cooling rates with the canonical dust-to-water ratio (see text), a = 0.0, 

T = 1500 K 



Nh,o 










•^LTEapprox 




% error 


1013 


2, 


.704 X 10' 


12 


2, 


.705 X 10' 


12 


0, 


.02 


1014 


2, 


.699 X 10' 


-12 


2, 


.702 X 10' 


12 


0, 


.12 


1015 


2, 


.664 X 10' 


-12 


2, 


.677 X 10' 


12 


0, 


.46 


1016 


2, 


.502 X 10' 


-12 


2, 


.527 X 10' 


12 


1, 


.01 


lOi^ 


2, 


.031 X 10' 


-12 


2, 


.055 X 10' 


12 


1, 


.36 


lOis 


9, 


.849 X 10' 


-13 


9, 


.430 X 10' 


13 


4, 


.26 


1019 


7, 


.223 X 10" 


-13 


6, 


.338 X 10" 


13 


12, 


.26 


1020 


3, 


.081 X 10" 


-17 


2, 


.612 X 10' 


17 


15, 


.23 


1021 


2, 


.210 X 10" 


-28 


4, 


.353 X 10' 


-28 


97, 


.01 


1022 




0.00 






0.00 




0, 


.00 


1023 




0.00 






0.00 




0, 


.00 


1024 




0.00 






0.00 




0, 


.00 
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Table 7. Cooling rates without the inclusion of dust, a = 0.0, T = 1750 K 



NhzO -^LTEexact -^LTEapprox ^ ^1101 



1 ni3 

iU 


A 

4. 


91 Q 

,ziy 


X 


iU 


-12 


A 

4, 


.zzu 


X 


lU 


-12 


n 

u, 


.Ui 




4, 


,214 


X 


10" 


-12 


4, 


.217 


X 


10- 


12 


0, 


.08 


10^5 


4, 


,177 


X 


10" 


-12 


4, 


.193 


X 


10- 


-12 


0, 


.37 




4. 


,002 


X 


10" 


-12 


4, 


.040 


X 


10- 


-12 


0, 


.94 


10^^ 


3. 


,501 


X 


10^ 


-12 


3, 


.575 


X 


10- 


-12 


2, 


.11 


10^8 


2. 


,446 


X 


10^ 


-12 


2, 


.511 


X 


10- 


-12 


2, 


.68 


10^9 


1. 


,158 


X 


10^ 


-12 


1, 


.189 


X 


10- 


-12 


2, 


.74 


1020 


4, 


,219 


X 


10- 


-13 


4, 


.343 


X 


10- 


-13 


2, 


.93 


1021 


1. 


,240 


X 


10- 


-13 


1, 


.279 


X 


10- 


13 


3, 


.15 


1022 


2. 


,863 


X 


10- 


-14 


2, 


.958 


X 


10- 


-14 


3, 


.35 


1023 


4. 


,913 


X 


10- 


-15 


5, 


.101 


X 


10- 


-15 


3, 


.83 


1024 


6. 


,215 


X 


10- 


-16 


1, 


.156 


X 


10- 


16 


5, 


.75 



Table 8. Cooling rates for the canonical dust-to- water ratio (see text), a = 0.0, 

T = 1750 K 



Nh20 -^LTEexact '^LTE.pp.ox % error 



1013 


4. 


,219 X 10^ 


-12 


4. 


,220 X 10^ 


-12 


0, 


.01 


1014 


4. 


,213 X 10^ 


-12 


4. 


,217 X 10' 


-12 


0, 


.08 


1015 


4. 


,173 X 10^ 


-12 


4, 


,188 X 10' 


-12 


0, 


.35 


1016 


3. 


,973 X 10^ 


-12 


4, 


,006 X 10' 


-12 


0, 


.82 


lOi^ 


3. 


,329 X 10^ 


-12 


3, 


,365 X 10' 


-12 


1, 


.15 


lOis 


1. 


,754 X 10^ 


-12 


1, 


,696 X 10' 


-12 


3, 


.27 


1019 


1. 


,383 X 10" 


-12 


1, 


,233 X 10' 


-12 


10, 


.84 


1020 


5, 


,735 X 10" 


-17 


4, 


,899 X 10' 


-17 


14, 


.57 


1021 


3. 


,973 X 10" 


-28 


7. 


,831 X 10' 


-28 


97, 


.13 


1022 




0.00 






0.00 




0, 


.00 


1023 




0.00 






0.00 




0, 


.00 


1024 




0.00 






0.00 




0, 


.00 
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Table 9. Cooling rates without the inclusion of dust, a = 0.0, T = 2000 K 



Nh20 -^LTEexact -^LTEapprox ^ ^1101 



1 ni3 

iU 


0, 


.Uoo 


X 


iU 


-12 


U- 


ns/i 

.Uo4 


X 


iU 


-12 


n 


.Ui 




6, 


.078 


X 


10" 


-12 


6, 


.082 


X 


10- 


12 


0, 


.02 


lO^s 


6, 


.040 


X 


10" 


-12 


6, 


.057 


X 


10- 


-12 


0, 


.28 


1016 


5, 


.847 


X 


10" 


-12 


5, 


.893 


X 


10- 


-12 


0, 


.78 


lOi^ 


5, 


.250 


X 


10^ 


-12 


5, 


.344 


X 


10- 


-12 


1, 


.79 


1018 


3 


.943 


X 


10^ 


-12 


4, 


.045 


X 


10- 


-12 


2, 


.58 


1019 


2, 


.084 


X 


10^ 


-12 


2, 


.139 


X 


10- 


-12 


2, 


.68 


1020 


8, 


.297 


X 


10- 


-13 


8, 


.536 


X 


10- 


-13 


2, 


.89 


1021 


2, 


.617 


X 


10~ 


13 


2, 


.698 


X 


10- 


13 


3, 


.12 


1022 


6, 


.260 


X 


10- 


-14 


6, 


.464 


X 


10- 


-14 


3, 


.26 


1023 


1, 


.013 


X 


10- 


-14 


1, 


.724 


X 


10- 


-14 


3, 


.83 


1024 


1, 


.060 


X 


10- 


-14 


L 


.118 


X 


10- 


-14 


5, 


.54 



Table 10. Cooling rates for the canonical dust-to-water ratio (see text), a = 0.0, 

T = 2000 K 



Nh,o 










•^LTEapprox 




% error 


1013 


6, 


.083 X 10- 


12 


6, 


.084 X 10- 


-12 


0, 


.01 


1014 


6, 


.077 X 10- 


-12 


6, 


.081 X 10- 


12 


0, 


.06 


1015 


6, 


.034 X 10- 


-12 


6, 


.051 X 10- 


12 


0, 


.28 


1016 


5, 


.803 X 10- 


-12 


5, 


.842 X 10- 


12 


0, 


.69 


lOi^ 


4, 


.977 X 10- 


-12 


5, 


.025 X 10" 


12 


0, 


.98 


1018 


2, 


.783 X 10- 


-12 


2, 


.716 X 10" 


12 


2, 


.39 


1019 


2, 


.351 X 10- 


-12 


2, 


.129 X 10- 


12 


9, 


.45 


1020 


9, 


.721 X 10- 


-17 


8, 


.359 X 10- 


17 


14, 


.01 


1021 


6 


.497 X 10- 


-28 


1, 


.283 X 10- 


-27 


97, 


.59 


1022 




0.00 






0.00 




0, 


.00 


1023 




0.00 






0.00 




0, 


.00 


1024 




0.00 






0.00 




0, 


.00 
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Table 11. Cooling rates without the inclusion of dust, a = 0.0, T = 2250 K 



Nh20 -^LTEexact -^LTEapprox ^ ^1101 



1 ni3 

iU 


S 9S/I 
o.Zo4 


X 


1 n- 

iU 


-12 


o.ZoO 


X 


iU 


-12 


n 


.Ui 




8.279 


X 


10" 


-12 


8.283 


X 


10- 


-12 


0, 


.04 




8.241 


X 


10" 


-12 


8.260 


X 


10- 


-12 


0, 


.22 




8.040 


X 


10" 


-12 


8.092 


X 


10- 


-12 


0, 


.65 


10^^ 


7.364 


X 


10^ 


-12 


7.476 


X 


10- 


-12 


1, 


.52 




5.835 


X 


10^ 


-12 


5.979 


X 


10- 


-12 


2, 


.47 




3.383 


X 


10^ 


-12 


3.471 


X 


10- 


-12 


2, 


.63 


1020 


1.457 


X 


10- 


-12 


1.499 


X 


10- 


-12 


2, 


.85 


1021 


4.870 


X 


10- 


-13 


5.022 


X 


10- 


-13 


3, 


.11 


1022 


1.185 


X 


10- 


-13 


1.222 


X 


10- 


-13 


3 


.20 


1023 


1.726 


X 


10- 


-14 


1.793 


X 


10- 


-14 


3, 


.87 


1024 


1.440 


X 


10- 


-15 


1.543 


X 


10- 


-15 


7, 


.14 



Table 12. Cooling rates for the canonical dust-to-water ratio (see text), a = 0.0, 

T = 2250 K 



Nh,o 








•^LTEapprox 




% error 


1013 


8, 


.284 X 10" 


12 


8.285 X 10" 


12 


0.01 


1014 


8, 


.279 X 10" 


-12 


8.281 X 10" 


12 


0.04 


1015 


8, 


.233 X 10" 


-12 


8.250 X 10" 


12 


0.22 


1016 


7, 


.976 X 10" 


-12 


8.022 X 10" 


12 


0.58 


lOi^ 


6, 


.963 X 10" 


-12 


7.023 X 10" 


12 


0.86 


lOis 


4, 


.065 X 10" 


-12 


3.998 X 10" 


12 


1.64 


1019 


3, 


.658 X 10" 


-12 


3.361 X 10" 


12 


8.11 


1020 


1, 


.532 X 10" 


-16 


1.324 X 10" 


16 


13.56 


1021 


9, 


.895 X 10" 


-28 


1.954 X 10" 


■27 


97.46 


1022 




0.00 




0.00 




0.00 


1023 




0.00 




0.00 




0.00 


1024 




0.00 




0.00 




0.00 
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is not clear is the degree to which hne coohng is significant in the region immediately past 
the shock. Does line cooling lead to a significant drop in temperature (> 10^ K) before 
dust grains begin to reabsorb line photons? Are the cooling rates of chondrules dominated 
by line cooling and therefore high (~ lO'^Khr"^) at the temperatures at which chondrules 
crystallize, as found by INSN, or can line cooling be neglected, as DC02 and CH02 assume? 
To answer these questions we have constructed a toy model to assess the maximum possible 
importance of line cooling and to determine whether more detailed calculations of chondrule 
formation in nebular shocks need to include line cooling. 

Our toy model builds on the calculation presented in Section 1. We find it convenient 
to convert a time derivative d/dt to a spatial derivative V^d/dz, and then convert the spatial 
variable z (distance past the shock front) into a column density of water past the shock front 
(assuming the water density remains constant). Then the cooling with column density of 
water is given by 




The cooling rate is a function of the water column density but also the dust-to-water ratio 
(via £lte)- Given a dust-to-water ratio, this equation can then be integrated to find T as 
a function of Wh20- Note that we have not accounted for the fact that the line broadening 
should decrease as the temperature drops, making it slightly harder for the gas to cool; our 
simplified analysis therefore overestimates the cooling somewhat. 

We plot the results of this integration in Figure [9] for both the dust-free case and 
the case where the dust is one-tenth that of the canonical value (in order to show the 
effect of even a small amount of dust). For ease of comparison to chondrules we have also 
assumed the canonical water-to-gas ratio (and other canonical values) to convert A^h20 to 
a time. Specifically, we assumed a pre-shock density of Pg = 10~^gcm~^, a shock velocity 
oi Vs = 7kms~^, assumed the density increased and the velocity decreased by a factor of 
(7 + l)/(7 — 1) = 6 past the shock front, and used our canonical water-to-gas ratio of 
'nH2o/nH2 = 8 X 10~^. These parameters are consistent with our assumption of an initial 
post-shock temperature of 2200 K (DC02; CH02). In the absence of dust, the gas would 
cool below 1400 K in roughly 300 seconds, leading to a cooling rate of > 1 x lO'^Khr"^. 
This cooling is attributable solely to effective cooling by water; by comparison to NK93 
we infer these are mostly vibrational photons. As chondrules and gas are expected to be 
thermally well coupled more than 100 seconds past the shock front (DC02; CH02), this can 
be interpreted as a likely cooling rate of chondrules as well. The cooling rate begins to taper 
off as large column densities are reached, but not until times > 10^ s, by which time the gas 
and chondrules have traveled roughly 10^'^ cm past the shock front, equivalent to a column 
density of water > 10^^ cm~^. 
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Fig. 9. — Gas temperature as a function of time, both with and without the inclusion of 
dust grains that can absorb hne photons (see text for details). 
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Including the absorption of line photons by dust grains completely changes the thermal 
history of the gas and chondrules. Using an opacity of dust ~ 30 cm^ g~^ (at short wave- 
lengths) and a post-shock gas density of 6 x 10~^gcm~^, the optical depth of dust exceeds 
unity after only 55 km past the shock front, which is reached after traveling only 48 seconds. 
During this first minute or so past the shock front, the line photons emitted by H2O molecules 
escape freely and effectively cool the gas; but after the first minute or so, they are absorbed 
by dust grains instead. Instead of H2O line photons freely escaping to the shock front and 
cooling the gas and dust system, the photons are reabsorbed "on-the-spot" and do not cool 
the gas. This effect of the reduced cooling rate is clearly seen in Figure [9] beyond about 10^ 
seconds. The trapping of line photons by dust grains becomes so effective, no further cooling 
by line emission is possible after a few minutes. Ignoring the effects of cooling other than 
line photons, the temperature would stabilize at about 1600 K. For this particular choice 
of parameters, chondrules would not cool at rates ~ lO^Khr"^ through their crystallization 
temperatures (1400 - 1800 K) as in INSN and Miura & Nakamoto (2006); in fact, they 
wouldn't cool at all under the assumptions of our toy model. 

We have investigated the effect of varying the dust-to-gas and water-to-gas ratios on 
the cooling rates of chondrules. Their behaviors are similar to those depicted in Figure M 
In Table 4 we report only the average cooling rates of chondrules over their crystallization 
temperature range (i.e., from 1800 K to 1400 K). The gas density and shock velocity are 
not varied. The results of this parameter study show that even a small amount of dust very 
effectively shuts down the line cooling due to H2O. 

4. Discussion and Conclusions 

We have built upon the line cooling work of NK93, by using a much expanded database of 
spectral lines, improved escape probabilities, and the inclusion of dust grains. We have found 
that although the cooling rate due to H2O line cooling is a complicated combination of gas 
and dust opacity, the cooling rate can be well approximated by i2LTE(T, Td) ~ '^^LTE(Td = 0) x 
£"2(0.95 Td). This allows us to investigate the maximum cooling effect H2O line cooling will 
have in chondrule-forming nebular shocks. 

Because vibrational line photons from H2O so effectively cool the gas, the high oxy- 
gen fugacities and H2O densities inferred for the chondrule-forming environment (Fedkin & 
Grossman 2007) would seem incompatible with cooling rates < 100 K/hr, but we find that 
dust significantly reduces the photon escape probabilities and the cooling rates. Dust grains 
absorb the line photons that would escape and cool the gas. Without dust, H2O line cooling 
reduces the temperature of the gas to 1400 - 1800 K in < 0.1 hr and cooling continues at this 
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Table 13. Cooling Rates for various dust to gas and water to gas ratios 



dust / gas^ 


water /gas^ 


T at 10^ s (K) 


Cooling Rate ^ 


0.0 


0.1 


1275 


2.7 X 10^ 


0.0 


1.0 


1000 


3.2 X 10^ 


0.0 


10.0 


875 


3.7 X 10^ 


0.1 


0.1 


2000 


N/A 


0.1 


1.0 


1650 


d 


0.1 


10.0 


1225 


2.3 X 10^ 


1.0 


0.1 


2175 


N/A 


1.0 


1.0 


2000 


N/A 


1.0 


10.0 


1600 


c 


10.0 


0.1 


2200 


N/A 


10.0 


1.0 


2175 


N/A 


10.0 


10.0 


2000 


N/A 



'^times the canonical dust/gas mass ratio 5 x 10 ^ 

Himes the canonical water/gas ratio rimo/nm = 8 x 
10-^ 

■^the average cooling rate between 1800 and 1400 K, in 
Khr"^ 

'^490Khr"^ between 1800 and 1650 K; never cooled 
below 1650 K 

^5000 Khr"^ between 1800 and 1600 K; never cooled 
below 1600 K 
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rate (>10^ K/hr). With dust, however, the grains absorb hne photons and inhibit coohng, 
leading to a cessation of coohng after ~ 0.1 hr (~ 8 min). Line coohng is a very effective 
and important coohng agent in the first few minutes and must be included in any compre- 
hensive shock model. After the first few minutes, however, the thermally well-coupled gas 
and chondrules cool together at a slow rate (only as fast as they travel ~ 1 optical depth 
through the dust), as in DC02 and CH02. INSN failed to note this effect; although they 
treat the gas and dust as well coupled (Tgas = Tjust), their system was effectively opticahy 
thin to line emission, allowing all line photons to escape the region and cool the gas. 

As we are interested at this time only in the maximum cooling effects of line photons 
(the "worst case scenario"), our calculations of cooling rates made no assumptions about 
radiative transfer; the chondrule temperatures followed the gas temperature. Thermal his- 
tories of chondrules including radiative transfer, line cooling, and other relevant effects will 
be reported on in future work (Morris et al. 2009b, in preparation). Since we have assumed 
that line photons either escape completely or are absorbed immediately, the true situation 
will clearly be bracketed between the two extremes shown in Figure [91 In these calculations, 
we have assumed that the dust is thermodynamically coupled to the gas. All temperatures 
(gas, dust, and chondrules), in actuahty, will be determined by dust /chondrule opacity and 
radiative transfer. DC02 showed that the cooling rates of chondrules (without the inclusion 
of line cooling) were given by the formula that appears above as Equation 1. In the case 
considered in Figure [9], the chondrule concentration may be considered to be low (C ^ 1), 
and the dust opacity is characterized hj 6 = 0.1 x 28.76/1.14 = 2.5. This yields a total 
cooling rate of gas and dust through the crystallization range of about 125Khr~^ after line 
cooling shuts off. Higher chondrule concentrations (C > 10^) appear to have been typical 
during chondrule formation (Cuzzi & Alexander 2006; Alexander et al. 2008). These would 
increase the cooling rate to several xlO^Khr"^. Line cooling is significant during the first 
few minutes past the shock, and more detailed calculations will be necessary to constrain 
the initial drop in temperature from the peak to below the liquidus. However, our results 
have shown that dust effectively shuts down line cooling within minutes and it then takes 
hours for chondrules to crystallize in a shock, consistent with DC02 and CH02. 

As mentioned previously, we did not calculate line cooling due to CO, as it is a factor 
of ~ 20 less effective than H2O. CO also cools by emitting line photons, but dust will 
absorb these line photons as readily as those emitted by H2O. Figure [9] shows a drop in 
temperature of ~ 600 K. If CO were included, the temperature would probably drop an 
extra 5% (~ 630 K), which is within the uncertainties of the H2O abundances we have 
assumed. 

This study has shown that the gas and chondrules will cool much too rapidly to match 
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the experimental constraints on chondrule cooling rates if dust is not present. The question 
of the dust abundance and its opacity therefore becomes paramount. But, significantly, we 
have found that enhancing only the amount of H2O causes higher cooling rates, whereas 
enhancing the amount of H2O and dust together does not. What affects the cooling rate 
the most, therefore, is not the amount of water, but the dust-to-water ratio. Potential 
complications, though, include the possibility that dust will evaporate immediately upon 
passage of the shock, and the possibility that dust vapor may recondense in the post-shock 
region (e.g., Scott & Krot 2005). 

In a related vein, we note an application of this work only tangentially related to chon- 
drule formation. Ciesla et al. (2003) suggested that phyllosilicates found in the fine-grained 
accretionary rims of chondrules in CM meteorites may have resulted from rapid gas-phase 
reactions of water vapor with silicates. Previous studies indicated that phyllosilicate pro- 
duction would be kinetically inhibited in the solar nebula (Prinn & Fegley 1987); however, 
Ciesla et al. (2003) showed that nebular shocks could result in a local increase in the water 
vapor pressure, thereby increasing the rate of phyllosilicate formation, as well as increasing 
the temperature at which phyllosilicates become stable. Chondrule- forming shocks in icy 
regions of the solar nebula could therefore account for both the formation of chondrules and 
their fine-grained phyllosilicate rims. A possible objection to this hypothesis is that the high 
abundance of water vapor in the shocked region would cool the gas too rapidly to allow 
the production of much phyllosilicates. As our work here shows, however, rapid cooling 
by line photons will not be associated with shocked, water-rich nebula gas, as long as the 
dust-to-water ratio is close to the canonical value for the solar nebula. An investigation of 
this hypothesis is also planned for future work. 

This work was supported by NASA Origins of Solar Systems grant NNG06GI65G. We 
would like to thank Frank Timmes for assistance on certain computing aspects of the prob- 
lem. 
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A. Line Cooling and Optical Depth 

The cooling per line is given by the expression 

Kul n. 



nH20 ^H20 



This preprint was prepared with the AAS IM^rpC macros v5.2. 
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where (in LTE) 

g-WfeT^ (A2) 



_ ni Qu Hv/kT 



nUiO nUiO 9i 



and the Einstein coefficient, A^i is given by 



^^J^^, (A3) 



Substituting Eq. IA2I and Eq. IA3I into Eq. lAll gives 



^ ^ e-W^^^ ( ^ B,^ 1 hu 



n-HiO \nH20 / V c- 
= -;— ^lu — e 

rii hu ^ , ^hu/kr\ ( 8nhu^ ^-hu/kT 



rii hviu 



R h .-/WfcTx /^Stt/^z/^ e-W^^ \ 
5^u {l-e-'-'/'^){4nB,{T)) (A4) 



What we are given in the SCAN-H2O hne hst is hnestrength. The 1.2 milhon hues 
of spectral data included in the line list cover a range in wavelength from ~ 6700 A to 
25 /xm. The data consist of the following: the wavenumber at line center, uq, in cm~^, the 
temperature-independent line strength, 5*0, in km/mol, and the excitation energy, Eiow, in 
cm~^. From the temperature-independent line strength, Sq, given in the list, the temperature- 
dependent hnestrength, S{T), is calculated as follows: 



g^j.^ ^ exp ( -E,,^-^) (1 - exp (-z/p^)) ^ ^^^^ 

Qvib{,T)Qrot{T) 



The vibrational partition function, Qmbi is given for temperatures ranging from 200 K 
to 8000 K, in increments of 200 K, in the documentation for the SCAN-H2O line list. We 
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have interpolated values for Qyn, using a quadratic spline. The rotational partition function, 
Qrot, is calculated with a subroutine provided in the documentation for the SCAN-H2O line 
list. S{T), as shown by Eq. IA5[ is in units of km/mol, and we wish to have S{T) in cm^ ■ Hz, 
so we must convert wavenumber to frequency which gives 



S{T,u) = 10' . So^, (A6) 



where A^^ is Avogadro's number, and we have converted from km to cm. This results in a 
conversion factor where 



S{T,u) =4.98 xlO~^ So = S{T). (AT) 
The left side of Eq. IA4I is just S{T), as given by Eq. IA51 so in terms of linestrength 



Aul ^^^^ ( Snhu' 



According to Plume et al. (2004), the optical depth to line center is given by 

To = A^i —- , (A9) 

1 /2 

where At; = (2kT/m) ' and the Einstein coefficients, A^i and Bi^ are given as follows: 

A„, = H^ff B„, = H^^:! £i (AlO) 



which gives the gas optical depth as 
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= A., ^ 1 4 ^ (1 ■ (All) 



In local thermodynamic equilibrium (LTE), 



iV. = iV^.o ^'!"P/Tg-/^^\ (A12) 



which gives 



Integrating over all frequencies, and normalizing to Ai^o = Af/c z/q gives the temperature- 
dependent hne strength 

5., = 5(T) = ^ ^ 5,. (1 - e-'^^/^^) , (AM) 
where the units are cm^ Hz molecule"^, resulting in the gas optical depth to line center 

To = Nn.o ri^KI\ ■ 



We want the frequency-integrated optical depth, however, which is given by 

(At;/c) I/O 



